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PREFACE 


In  pulse-jets  of  the  type  used  to  propel  the  German 
Buzz  Bomb,  turbulence  appears  to  have  considerable  influence 
on  combustion  and  on  effective  flame  speeds  and  therefore 
on  the  oscillation  and  flow  of  the  gas  through  the  motor. 
The  following  report  assumes  that  these  effects  can  be 
accounted  for  mathematically  by  use  of  parameters  (to  be 
determined  by  controlled  fundamental  experiments)  in 
general  one-dimensionalized  thermodynamic al  gas  equations. 
Idealized  numerical  cases  are  considered.   It  is  intended 
that  these  equations  be  used  to  investigate  optimum  internal 
conditions  which,  it  is  hoped,  will  lead  to  specific 
recommendations  for  the  design  of  more  effective  pulse- 
jets,  as  well  as  to  a  fundamental  understanding  of  auto- 
matic reignition  and  other  characteristic  phenomena. 
Previous  attempts  at  a  theory  of  intermittent  jet  motors 
have  ignored  flame  or  wave  propagation  as  well  as 
turbulence  effects. 

Professor  J.  K.  L.  MacDonald,  the  author  of  the 
present  report,  was  greatly  assisted  by  Dr.  S.  A.  Schaaf 
in  particular  and  by  other  members  of  the  staff  of  the 
Institute  for  Mathematics  and  Mechanics  in  Nev/  York  University. 
This  important  investigation  has  been  carried  out  in 
connection  with  a  research  program  under  Contract  No. 
N0a(s)-7370  with  the  Bureau  of  Aeronautics  of  the  Navy 
Department. 

R.  Courant 

Director 

Institute  for  Mathematics  and  Mechanics 
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k  GAS  DYNAMICAL  FORMULATION  FOR 
WAVES   AND  COMBUSTION  IN  PULSE- JETS 


SUMMARY 

A  brief  outline  Is  given  of  known  experimental  data 
on  the  motions  of  the  gases  and  flames  in  pulse-Jets  (such 
as  the  German  "Bussz-bomb"  motor).  The  considerable 
assumptions  made  in  previous  theoretical  treatments  are 
Indicated  and  a  rather  comprehensive  one-dlmensionallzed 
hydro-thermodynamlcal  formulation  for  pulse-Jets  is  then 
developed.  Terms  are  Included  corresponding  In  particular 
to:   (l)  wave  and  flame  propagation  throughout  the  jet, 
(11)  effective  duration  of  burning  in  the  turbulent  gases, 
and  (ill)  non-uniform  cross-section  and  boundary  conditions 
at  the  intake  and  at  the  exhaust  ends.   An  outline  Is  given 
of  the  nature  of  further  experimental  data  required  for  the 
explicit  determination  of  boundary  conditions  and  of  the 
parameters  Involved  in  combustion  and  flame  propagation  in 
a  tTirbulent  mixture.   It  is  strongly  xirged  that  such  data 
be  obtained  In  order  to  make  detailed  applications  of  the 
theory  to  practical  questions  of  efficiency  and  design.   An 
Idealized  example  is  solved  explicitly  (analytically  to  first 
order  terms)  on  the  basis  of  presumably  reasonable  assumptions 
concerning  such  data,  and  is  seen  to  represent  observed 
featvires  and  to  contradict  some  assumptions  made  in  previous 
theoretical  treatments.  In  particular  pressure  effects  due 
to  the  moving  fleme  fronts  appear  to  be  Important. 

A  variety  of  alternative  forms  is  presented  to  provide 
a  framework  into  which  future  observations  and  explicit 
examples  may  fit.  It  is  hoped  that  the  existence  of  such 
a  framework  will  help  to  stimulate  really  fundamental 
experimental  rese«u?ch  on  pulse-Jets. 
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A  GAS  DYNAMICAL  FORMULATION  FOR 
WAVES  AND  COMBUSTION  IN  PULSE- JETS 


1.   Introduction 

Following  the  development  of  the  V.l  "Buzz  Bomb"  by 
the  Germans,  various  groups  In  England  and  the  U.S.A.  started 
general  Investigations  of  pulse-Jet  motors.   Among  these  groups 
were  the  following:  The  R.A.E.  (Royal  Aircraft  Establishment, 
parnborougih,  England),'  '   the  U.S.  Army  Air  Force  (at  Wright 
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(a)  German  V.l  type  motor  with  axial  symmetry.  Frequency  about  50  cycles/se 
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(b)  Annapolis  "Resojet"  type  motor.  Frequency  about  100  cycles/sec. 


Pig.  1:  Approximate  structure  and  dimensions  of 
standard  pulse-Jets. 
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Pield),^    '   the  N.A.C.A.    (National  Advisory  Council  for 

(2)  (3) 

Aeronautics,  Cleveland),^  '  a  BuAer  group  (at  Annapolis),^  ' 

(4) 
the  Aerojet  Corp.  in  California,^  '  the  Reaction  Motors 

Corp.  in  New  Jersey,  the  Jet  Propelled  Missiles  Panel  of 

the  O.S.R.D. ^^^  (for  the  U.S.  Navy)  and  the  AMG-NYU  (Applied 

Mathematics  Group  at  New  York  University). 

The  operation  of  standard  types  of  pulse  jets,  such 

as  are  represented  in  Figure  1,  is  roughly  as  follows: 

The  motor  is  placed  in  an  air  stream  which  provides  a  "ram 

pressure"  to  force  air  through  the  reed  valves  (see  Pig.  1) 

into  a  region  where  the  air  is  mixed  with  a  continually 

sprayed  fuel  (such  as  gasoline).  The  mixture  is  initially 

fired  by  means  of  a  spark-plug,  and  the  btirning  raises  the 

gas  pressure  in  the  combustion  chamber,  ultimately  causing 

the  reed  valves  to  close  and  exhaust  gases  to  flow  out 

through  the  tail  pipe,  a  thrust  being  developed  during  the 

process.  The  exhaust  generates  rarefaction  waves  which 

soon  reduce  the  pressure  behind  the  reed  valves  until  they 

open  again  and  admit  more  air.  Then,  even  without  further 

use  of  the  spark-plug  a  self  firing  occurs  in  the  gas, 

apparently  in  a  central  region  near  the  section  through  the 

spark  plug,  and  a  regular  cycle  is  established.  The  mechanism 

of  the  self -firing  is  not  clearly  understood  but  it  seems 

probable  that  pressure  waves  reflected  from  the  ends  of  the 

motor  cause  a  local  heating  and  mixing  of  the  gas  sufficient 

for  its  ignition  near  the  interface  between  the  incoming 

(unburnt)  mixture  and  the  previously  burnt  gases.  The  flame 

apparently  spreads  from  the  region  of  re-ignition  with  a 

relatively  slow  speed,  in  a  manner  indicated  in  Figure  5; 

the  process  thus  resembles  a  wave  propagation  rather  .than 

an  instantaneous  explosion,   (In  an  ordinary  internal 

combustion  engine  the  distances  involved  in  wave  and  flame 

propagation  are  relatively  small  and  the  state  of  the  gas  is 
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effectlvely  homogeneous  at  all  times. ) 

Using  fairly  crude  instruments  most  of  the  above 
mentioned  groups  of  investigators  have  made  measurements 
of  seme  sorts  of  average  values  (over  the  cycle)  of  the 
following  quantities:  thrust,  drag,  ram  pressure,  rates 
of  flow  of  fuel  and  air,  and  pressures  and  temperatures  at 
various  points.   Some  attention  was  devoted  to  the  effects 
produced  by  changes  in  structures  and  in  fuels.   A  few 
nominal  pressure-time  oscillograms  were  obtained  but  were 
subject  to  unknown  perturbations  due  to  violent  thermal, 
mechanical  and  turbulent  oscillations  in  and  around  the 
pick-ups.   Preliminary  studies  of  the  motions  of  flames 
and  of  exhaust  gases  were  made  at  Wright  Field,  at  Cleveland 
and  at  Annapolis,  by  use  of  high  speed  movie  cameras 
(2000  frames  per  second).  Really  fundamental  and  accurate 
experimental  work  remains  to  be  done.  Under  contract  to 
BuAer,  groups  at  the  N.H.L.  (Naval  Research  Laboratory, 
Washington),  at  the  Aerojet  Corporation  and  at  New  York 
University  hope  to  carry  through  fairly  comprehensive 
programs  of  experimental  work  on  the  "internal  ballistics" 
of  pulse-jets.  The  results  of  such  investigations  would 
provide  basic  data  for  use  in  a  theoretical  treatment  such 
as  that  which  will  be  discussed  in  section  3,  and  which, 
it  is  hoped,  will  lead  to  improvements  in  the  designs  of 
pulse-jets.  The  nature  of  desired  experimental  data  is 
indicated  in  Section  4.   It  is  strongly  urged  that  such  data 
be  obtained  by  groups  with  the  required  facilities. 

2.   Previous  Theoretical  Treatments  of  Waves  in  Pulse-Jets. 

Prom  the  few  relevant  reports  received  by  the 
AMG-NYU  it  would  anpear  that  British  observers  first  noted 
that  the  frequency  of  the  V.l  motor  under  normal  operating 
conditions  is  almost  the  same  as  that  of  the  fundamental 
mode  in  an  organ  pipe  having  the  same  dimensions  and  the  same 


-4- 

( average)  gas  temperature  as  the  motor. 

In  a  Wright  Field  report^  the  V.l  motor  is  treated  as 
a  sort  of  Helraholtz  resonator. ^   The  usual  Helmholtz 

resonator  has  roughly  the 
shape  of  a  flask,  as  indicated 
in  Pig.  2,  with  the  volvme  V 
and  its  cross-section  con- 
siderably greater  than  those 
Pig-   2  of  v.   Then  V  acts  like  a 

spring  against  which  oscillates 
the  mass  of  V*  (plus  small  effective  masses  at  the  ends  of 
V*)»   In  the  Buzz  Bomb  motor  V  and  V*  are  practically  un- 
distinguishable  from  one  another  and  the  asstunptions  used 
in  the  standard  resonator  theory  are  violated.   It  is  most 
probably  merely  a  numerical  accident  that  the  observed 
frequency  of  the  standard  Buzz  Bomb  is  in  close  agreement 
with  the  theoretical  resonator  frequency,  especially  since 
the  Wright  Field  report  treats  V  +  V  as  the  spring  for  V*. 
A  considerably  refined  (but  similar)  treatment  is  given  in 
two  reports  '  by  Professor  Nell  P.  Bailey  and  his  associates: 
The  pulse-jet  is  again  treated  mathematically  on  the 
assumption  that  the  state  of  the  gas  in  the  combustion 
chamber  remains  uniform  as  it  varies  with  time  and  produces 
forced  oscillations  of  acoustic  type  in  the  tail  pipe.   Some 
effects  due  to  non-instantaneous  burning  of  the  fuel  are 

included  by  assuming  a  periodic  exponential  pulse  form  for 

1  d© 
the  quantity  -g  gr  »  where  ©  is  the  absolute  temperature  in 

the  combustion  chamber  at  time  t.  Rather  comprehensive 

numerical  calculations  were  made  by  Bailey  but  the  assumption 

of  uniformity  of  state  (i.e.  no  wave  motion)  throughout  the 

combustion  chamber,  as  well  as  assumptions  about  the  bulk 

modulus  and  the  teraperatuie  of  the  gas,  seem,  on  the  basis 

of  a  more  complete  theory  to  have  little  Justification. 
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In  an  infonnal  meraorandxim^  by  Professor  H.  Lewy  of 

the  AMG-NYTT  there  is  indicated  an  interesting  type  of 
forimilation  for  waves  in  pulse-jets.   Although  Lewy's 
formulation  neglects  the  non-instantaneous  character  of  the 
combustion  and  makes  some  assumptions  of  doubtful  validity 
concerning  the  velocities  and  densities  In  the  gases,  it 
may  lead  to  useful  generalizations  after  more  experimental 
data  are  available.  The  method  depends  upon  the  effective 
elimination  from  the  (one  dimensional)  equations  of 
motion  and  continuity,  of  the  position  variable  x.  This 
elimination  results  from  an  x-integration  of  both  sides  of 
the  equations,  coupled  with  the  assumptions  that 

/  u(x,t)dx/u(L,t)  and  /  ^  (x,t  )dx/yo(0,t )  are  both 
'^o  ''o 

independent  of  the  time  t,  and  that  the  entropy  of  the  gas 

remains  constant.   Here  u(x,t)  is  the  gas  velocity  and 

/o(x,t)  the  gas  density,  and  x  =  0  and  x  =  L  respectively 

correspond  to  the  positions  of  the  intake  and  the  exhaust 

ends  of  the  pulse-Jet.   Non-linear  differential  equations 

then  arise  for  u(L,t)  and  /C(0,t)  and  are  solved  (for 

special  initial  conditions)  by  a  combination  of  analytical 

and  numerical  methods.   Generalizations  of  Lewy's  method, 

to  include  burning  periods  and  justified  relations  between 

averages,  are  planned  in  future  investigations  when  basic 

experimental  data  are  available. 

German  articles  covering  their  war  research  on 

theories  of  pulse-jets  are  now  becoming  available  (see 

reference  IS}.   All  such  articles  thus  far  received  have 

not  considered  flame  motion  and  non-instantaneous  combustion, 

but  they  have  given  extensive  strictly  nuirierical  analyses 

of  non-linear,  isentroplc  v/ave-motion  in  open  and  semi-open 

tubes  of  variable  cross-section. 
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3,      One-dimensionallzed  Formulation  for   Pulse-jets. 


Pi        r—, i '2 
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Pig.  3   Characteristic  shape  of  the  interior 
of  pulse-jets. 

p,  =  ram  pressure,  pg  =  atmospheric 

press\ire  out- 
side the  exhaust 
pipe. 


The  usual  pulse-jet  is  a  steel  tube  of  mildly  non-uniform 
cross-section  (circular  or  rectangular)  with  the  exhaust 
end  (see  Pig.  3)  constantly  open  and  with  the  size  of  the 
relatively  small  openings  at  the  intake  end  varying  with 
time  in  a  way  depending  on  the  varying  effective  pressure 
difference  "between  the  two  sides  of  the  intake,  and  also 
depending  on  the  constants  of  the  reed  valves. 

Although  detailed  observations  are  still  lacking  it 
seems  definite  that  in  the  pulse-jet  there  is  considerable 
turbulence,  probably  traceable  mainly  to  the  reed  valves 
and  to  inhomogeneity  in  the  fuel-air  mixture.   Correspondingly 
the  combustion  appears  to  proceed  inliom.ogeneously  in 
cluster  form  rather  than  in  the  relatively  slowly  moving 

idealized  plane  sheet  form  to  which  the  LeChatelier  theory 

(8 

and  its  recent  extensions  apply. ^   In  fact  it  is  not  unlikely 
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that  the  flsmes  spread  in  the  forra  of  turhulently  distorted 
and  interpenetrating  thin  shells  which  make  entire  regions 
of  the  gas  glow  with  conihustion  for  several  iriilliseconds 
(see  Pig.  5).   Presuinably  the  turbulence  greatly  increases 
the  effective  rate  of  diffusion  of  heat  and  of  chain- 
carriers  and  thereby  produces  the  observed  high  rates  of 
spread  of  the  flsjnes.  Locally  the  combustion  of  particles 
Is  probably  almost  instantaneous  after  their  ignition  by 
heat  and  chain-carriers  from  adjacent  particles,  but  in  a 
given  layer  of  the  swirling  mixture  of  thickness   A x  at  x 
(see  Pig.  3)  only  a  small  fraction  of  the  molecules  may 
undergo  combustion  during  any  small  interval  of  time   At» 
Correspondingly,  locally,  the  concepts  of  thermodynamic 
state  may  not  apply,  but  the  concepts  of  average  thermo- 
dynamic state  and  of  effective  non-ssero  duration  of  burning 
would  be  expected  to  have  a  meaning  for  the  layer  as  a 
whole.  The  term  "flame  speed"  could  be  applied  to  the 
effective  rate  of  advance  of  combustion  from  layer  to  layer 
along  the  x  direction.  Prom  this  viewpoint  it  would  seem 
quite  probable  that  the  cross-sectional  averages  of  salient 
physical  quantities  can  be  closely  approximated  by  use  of 
a  one-dimensionalized  model  (with  refinements  awaiting  the 
development  of  adequate  theoretical  treatments  of  turbulence 
and  combustion,  as  well  as  adequate  experimental  studies). 
In  such  a  one-dimensionalized  model  the  state  of  the  gas 
is  taken,  in  effect,  as  uniform  over  the  cross-section  S(x) 
(see  Fig.  3),  and  the  following  symbols  may  be  introduced 
(typical  consistent  units  are  indicated  for  use  in  subsequent 
equations):  At  distance  x(cm)  from  the  intake  at  time 
t(sec)  let 

a  =  Lagrange  coordinate,  i.e.  the  value  of  x  at  t  =  0  for 

the  particle  (layer)  which  is  at  x  at  time  t  (a  therefore 
varies  from  particle  to  particle). 
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m  =  'pseudo-Lagrange  coordinate*  with  the  following  meaning: 
For  the  particle  (layer)  which  reaches  the  position  x 
at  time  t,  ro  denotes  the  total  mass  of  gas  that  has 
flowed  thjpough  the  intake  from  some  convenient  reference 
tirne  t  until  the  time  of  arrival  t(m)  of  that  particular 
particle  at  the  (inside  end  of)  the  intake.   Each 
particle  therefore  has  its  own  value  for  m. 

u  =  velocity  (cm/sec)  of  the  gas  relative  to  the  tube. 

P  =  velocity  (cm/sec)  of  flsme  front  (at  x)  relative  to  the 
tube. 

p  =  pressure  (dynes/cm  static)  of  the  gas. 

yo  =  density  (gra/cm  )  of  the  gas. 

T=  l/o  =  specific  volume  (cm  /gm). 

©  =  absolute  temperature  (^Kelvin)  relative  to  the  moving  gas, 

ft  =  effective  rate  of  release  of  heat  energy  per  unit  mass 
of  gas  mixture  (ergs/gm  sec)  (at  x  at  time  t),  that  is 
the  sum  of  the  rate  of  release  of  chemical  energy  and 
of  the  rate  of  inflow  (per  unit  mass)  of  heat  by  conduction 
and  convection  from  adjacent  regions.  The  values  of  E 
should  be  determined  from  experiments  such  as  are 
indicated  in  section  4. 

TJ  =  internal  energy  (ergs/gm)  of  the  gas. 

C  and  C  =  specific  heats  (energies)  (ergs/gm *^K)  at  constant 
vol\ime  and  constant  pressure  respectively. 

Y  =  C  /C  =  ratio  of  specific  heats. 

R  =  C  -C  =  "ideal  gas  constant"  (ergs/gm°K). 
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3.1  Hydro-thermodynamlcal  equations  for  the  one-dinienslon- 
allzed  pulse-.jet. 

Note  that  the  above  definition  for  R  does  not  refer 

to  the  customary  'mole'  unit  but  is  inversely  proportional 

to  the  average  molecular  weight  of  the  mixture  at  the 

position  and  time  concerned.  For  commonly  used  fuel-air 

mixtures  R  remains  nearly  constant  but  y-l  drops  from  about 

.4  to  .2  as  e  rises  from  0°  to  about  2500*^0  during  the 

combustion.   Moreover,  for  such  mixtures  the  equation  of 

— ~ 

state  is  commonly  held^  to  be  represented  fairly  accurately 

by 

(1)  pT    =  R©  (R  constant), 
and  the  change  in  internal  energy  by 

(2)  du    =  c^de  =  ^lEfl 

In  ordinary  gas  dynamics  without  combustion,  the 
adlabatic  relation 

(3)  -pdr  =  C  d©  =  -^^r^   (whence  p  t''^  =  const,  for  y  const.) 

Is  usually  employed,  and  it  can  be  derived  from  the  first 
law  of  thermodynamics  and  the  relations  (1)  and  (2),   In 
the  present  formulation  the  symbol  E  is  introduced  to 
represent  all  heat  contributions.   (l),  (2)  and  the  first 
law  of  thermodynamics  then  yield  the  following  analogue  of 
the  adlabatic  relation  (3): 

(4)  Edt  -  pdr  =  C^d©  =  ^^^ 
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Here  the  differentials  refer  to  the  change  in  the  thermo- 
dynamic state  of  a  given  moving  particle  (layer)  during 
the  time  interval  dt.   For  the  case  !&  =  0  (4)  reduces  to 
r?).   An  enthalpy  equation  can  be  obtained  by  combining 

(4)  with  an  equation  of  motion  such  as  is  given  below,  or 
alternatively  by  considering  the  entire  flux  of  energy 
into  a  given  particle  (layer).   As  will  be  explained  in 
Section  3.2  the  introduction  of  shock  discontinuity 
conditions  is  apparently  an  unnecessary  refinement. 

There  are  various  pairs  of  independent  variables 
available  for  specifications  in  pulse-jets.  Eulerian 
coordinates  (x,t)^   are  probably  most  convenient  for 
discussing  steady  flows,  but  in  non-steady  flames  they 
lead,  in  general,  to  rather  cumbersome  relations. 
Lagrangian  coordinates  (a,t)^   keep  track  of  each  particle 
and  are  especially  convenient  for  treating  waves  and 
combustion  in  initially  resting  gases  in  the  tube;  but 
inflows  lead  to  some  complications  which  can  be  avoided  by 
use  of  pseudo-Lagrangian  coordinates  (m,t)  (v/here  m  is  as 
defined  on  page  8).  The  complications  concerned  are  due 
to  the  fact  that  the  Lagrange  coordinate  a  of  a  particle 
which  enters  the  intake  after  t  =  0,  depends  sensitively 
upon  the  dynamics  of  the  flow  outside  the  motor. 

In  Eulerian  coordinates  (x,t) 

(5)  S(x)/o^  (x,t)  +  Cs(x)/o(x,t)  u(x,t)]^  =  0 

is  the  equation  of  continuity,  and 

S(x)[yo(x,t^  u(x,t)]^  +  CS(x)  /o(x,t)  u^(x,t)]^  +  S(x)  p^(x,t)  =  0 
(6)-^  or,  alternatively 

^u^(x,t)  +  u(x,t)  u^(x,t)  +  p^(x,t)//o(x,t)  =  0 
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ls the  equation  of  motion.  Here  subscripts  x  and  t  represent 
corresponding  partial  derivatives  and  S(x)  is  as  defined 
In  Fig.  3  (page  6).   By  interpreting 

^  as  ^  +  u(x,t)^  , 

equation  (4)  can  be  expressed  In  terms  of  Eulerian 
coordinates.  Thus: 

(7)  (Y-l)E  =  YP(^>  u^)t  +  -rC^  4.  u^)p 

In  Lagrangian  coordinates  (a,t)  with  x(a,t)  =  a  at  t  =  0 
the  hydro-thermodynamlcal  equations  (4),  (5)  and  (6) 
become  respectively 

(8)  E  -  p(a,t)  r^(a,t)  =  C^  e^(a,t)  =  [p(a,t)  T(a,t)]^/(T-l), 

(9'  -at'.*)  =  l&lllll] 

and 

(10)  x^^  =-S(x)  T(a,0)  Pg^(a,t)/S(a). 

In  pseudo-Lagrangian  coordinates  (m,t)  the  corresponding 
equations  (see  Appendix  l)  are: 

E  -  p(m,t)  T^(m,t)  =  C^  ©^(m,t)  =  Cp(m,t)  T(m,t )]^/(y-1), 

(11)  whence 

(y-1)E  =  YPT^  *'^'Pt   ' 
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(12)  x^(m,t)  =  -r(in,t)/S(x) 

and 

(13)  x^^(Tn,t)  =  Six)   p^(ni,t) 

where  m  denotes  the  mass  of  the  mixture  that  has  flowed 
through  the  intake  from  some  reference  time,  t  ,  until  the 
time  t(m)  of  entry  there  of  the  particles  which  (in  the 
one  dlmensionalized  model)  will  reach  x  at  time  t. 
Symbolically 

t(m) 
m=/   [S(0+)  x^/T]dt  ^Q^ 

\ 

where  s(0+)  is  the  cross-section  just  following  the  valves 
at  X  =  0.   At  the  intake  end  the  mass  flux  is  given  by 

(14)  S(x)  x^(m,t)/T(m,t)  -»  -^^   for  x(m,t)  -»  0+  , 

with  j^»  ■■»  of  course  zero  during  intervals  when  the  reed 
valves  are  shut. 

3.2  On  approximations  for  use  in  solving  the  one-dimenslon- 
alized  equations. 

(i)   Shock  front  conditions.   It  is  well  known^   that 
compression  waves  (as  well  as  explosions)  can  give  rise  to 
shock  fronts  across  which  entropy  and  other  physical 
quantities  take  jumps  in  values.   It  might  therefore  be 
feared  that  a  treatment  of  pulse-jets  which  implicitly 
ignores  the  possibility  of  shock  waves  would  be  of  rather 
limited  value.   Such  a  fear  appears  to  be  groundless  when 
the  pressure  oscillations  and  the  burning  rates  ara  as  small 
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es  rough  measurements  in  standard  pulse-jets  Indicate.  Under 
such  circumstances  any  shocks  that  may  arise  in  the  gas 
mixture  would  he  expected  to  be  so  weak  that  p,u  and  t  in 
the  one-dimensionalized  model  could  be  regarded  as  continuous 
or  as  limits  of  continuous  approximations.   In  fact  if  one 
compares  pressure  distributions  calculated  respectively 
on  the  basis  of  (linear^  acoustic  wave  and  (non-linear) 
shock  wave  equations,  for  a  problem  such  as  that  indicated 
in  Figure  4,  one  finds  only  a  few  percent  difference  between 
the  two. 


Pi 


Pq  =  Pl/2 


p 

i 

Pi 

Po 
0 

RAREFACTION   SHOCK 
FRONT         FRONT 


V 


-^  X 


Time  t  =  0 


Time  t  >  0 


Fig.  4  ?/aves  produced  in  a  tube  in  which 
there  is  an  initial  discontinuity 
in  the  pressure  p  of  the  gas. 


The  omission  of  shock  conditions  permits  the  development 
of  relatively  simple  iteration  schemes  for  solving  the 
equations  analytically  rather  than  purel?/-  nvunerically,  a 
feattire  which  is  desirable  from  a  practical  viewpoint  because 
of  the  m.any  parameters  whose  influence  must  be  estimated 
for  desifn  and  other  purposes. 
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(ii)  Intake  and  exhaust  end  conditions*   At  the  exhaust 

end  (x  =  L,  see  Fig.  3),  rough  pressure  measurements 

fl2 

coupled  with  notions  used  in  acoustic  theory^   suggest  that 

to  a  good  approximation 

(15)  p  =  Po  at  X  =  L 

where  p^  is  the  effective  atmospheric  pressure  outside  the 

*  •  f  12 

exhaust  pipe  and  where  L  -  L  is  an  "open  end  correction"^ 

which  is  of  the  order  of  half  the  diameter  of  the  exhaust 

pipe.  With  this  ass-umption  the  effective  length  of  the 

pulse-Jet  is  L  . 

At  the  intake  or  reed-valve  end,  x  =  0,  suitable 

approximations  for  the  boundary  conditions  are  apparently 

not  so  obvious,  and  further  experimental  work  is  required 

to  nrovide  reliable  relations  (see  Section  4).   Since  the 

usual  reed  valve  openings  permit  flow  through  at  most  a 

small  fraction  of  the  cross-section  area  at  x  =  0  (see  Fig. 

3)  the  external  effective  ram  pressure,  p,,  probably  remains 

fairly  constant  during  a  cycle.   Presumably  when  the  internal 

pressure  p  at  x  =  0  exceeds  an  empirically  determinable 

value  p-  (depending  on  p,  and  on  the  constants  and  adjustment 

of  the  reed  valves)  the  valves  will  stay  shut  and  the  end 

condition  will  be 

(16)  u  =  0  when  p  >  p^  at  x  =  0+. 

(closed  valve  condition) 

V.Tien  p  <  p,  the  valve  openings  will  presumably  act  as  nozzles 

with  the  varying  flovr  through  them  satisfying  a  suitably 

( 13 
generalized  "pressure -he ad  loss"  relation^   ,  possibly  of  the 

form 
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(17)       h(p,u,p^)  =  0  for  p  <  Pj  at  x  =   0+. 

If  there  were  either  no  or  steady  combustion,  then  p  and  u 
in  (17)  would  have  steady  values,  say  p  and  u  .   If  the 
actual  oscillations  in  the  pulse- jet  are  regarded  as 
perturbations  of  the  steady  state  (as  will  be  explained  in 
(iv)  below)  then  at  x  =  0+ 


(18)         P  =  Pq  +  ^P  and  u  =  u  +  Au 


where  Ap  and  Au  are  the  perturbations  of  the  steady 
values  p  and  u  respectively.  For  sufficiently  small 
Ap  and  Au  (18)  in  (17)  yields  the  following  approximate 
boundary  condition: 

(19)        Ap  =  -z  Au  for  p  <  Pg  at  x  =  0+ 

(open  valve  condition)   with 

z  =  +[(9/au^)  h(p^,u^,p^)]/[(5>/8pQ)  h(Po,u^,p^)] 

a  constant  which  might  be  fairly  easily  determined  empirically 
if  (19)  is  an  appropriate  boundary  condition,   (Note  that 

(19)  is  equivalent  to  an  impedance  relation,  which  fact  would 
be  sufficient  justification  for  it  in  the  minds  of  many.) 

In  the  special  case  of  initially  steady  flow  it  would  be 
justified  to  take 

(20)  P  =  Pq  and  u  =  u   at  X  =  0+ 

(i.e.  just  inside  the  head  of  the  motor)  until  the  time  of 
arrival  of  compression  waves  from  the  initial  region  of 
ignition,  and,  asstiming  that  the  rise  in  p  at  that  time  is 
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so  rapid  as  to  be  considered  instantaneous,  thereafter  to 
take  (at  x  =  0+) 

(20' )    u  =  0   until  the  exhaust  waves  arrive  at  x  =  0* 


(ill)  Combustion  and  flame  speeds.  Figure  5,  based  on  an 
analysis   of  high  speed  movies  taken  through  holes  in  the 
side  wall  of  a  standard  V. 1  motor  (under  normal  conditions 
of  operation),  indicates  how,  during  a  single  cycle,  the 
visible  flames  spread  and  persist  after  self-ignition  at  the 
position  marked  I  in  Pig.  5. 


t  (sec) 


.02 


I 


NO  VISIBLE  FLAME 
DURING  THIS  TIME 


ONE  PERIOD   (about  .02  sec.) 


.01 


MINOR  BAND  OP 
'FLAME  PROBABLY  DUE 
?0  A  BRIEF  REBURNINO 
'OF  THE  ALMOST  COMPLETELY 
'BURNT  MIXTURE  LEFT  OVER 
''FROM  THE  PRECEDING  CYCLE 


10' 


-♦•  X 


Pig.  5   Flame  spread  during  a  typical 
cycle  of  the  V.l  motor.  The 
shaded  regions  denotd  burning. 
(The  gases  remain  hot  but 
normally  non-luminous  after 
combustion. ) 
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The  point  I  of  self -ignition  of  the  gas  probably  marks  the 
interface  between  the  incoming  fresh  mixture  and  the  almost 
completely  burnt  (but  chain-carrier  rich)  mixture  left  over 
from  the  proceeding  cycle.  Dud  to  the  motion  of  the  gases 
towards  the  exhaust  end  this  Interface  moves  with  the  particle 
speed  |ul ,  as  indicated  by  the  curve  II 'I",  while  the  minor 
band  IB  of  reburning  moves  with  a  flame  speed  | P|  (relative 
to  the  motor)  equal  to  the  s\ira  of  the  particle  speed  |u| 
and  the  'intrinsic  flame  speed'  f  relative  to  the  gas. 
On  the  basis  of  this  interpretation  a  numerical  analysis 
of  experimental  curves  corresponding  to  Pig.  5  indicates 
that  both  I  u|  and  f  amount  to  a  couple  of  hundred  feet  per 
second  at  I',  and  that  the  flame  speed  | P|  relative  to  the 
motor  can  be  as  large  as  2000  ft/sec  at  B.   (The  local  speed 
of  sound  at  B  is  probably  about  2500  ft/sec. )  The  maximum 
duration  of  combustion,  as  judged  from  luminosity,  does  not 
appear  to  exceed  .014  sec.  Doubtless  the  flame  pattern 
would  vary  with  the  conditions  of  operation  and  with  the 
geometry  of  the  tube.   In  fact  for  the  shorter  'Resojet' 
motor  represented  in  Pig.  1(b)  (page  1)  the  burning  has 
been  found  to  persist  for  a  larger  fraction  of  the  period 
(.01  sec),  and  to  be  not  so  sharply  outlined  as  that 
indicated  in  Pig.  5,   The  Resojet  is  therefore  probably  not 
very  suitable  for  fundamental  investigations  which  require 
consistently  stable  cycles. 

To  a  first  approximation  the  effective  rate  E  (see  p.  8) 
of  release  of  heat  energy  in  any  given  layer  of  the  mixture 
is  given  by 

E  =  0  before  and  after  combustion 
(21) 


=  AE/At  during  combustion. 
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where    At   is  the  duration  of  the  burning  and    AB  is  the 
energy  of  combustion  per  gm  of  the  mixture.      (For  a 
gasoline-air  mixture  in  the  weight  ratio  1:16,  with 
At  =  .014  sec,    AE/At  =  2.10       cm  /sec   . )     On  the  basis 
of  chemical  kinetics  E  might  be  expected  to  increase  with 
both  the  temperature  and  the  concentration  of  the 
combustible  part   of  the  mixture.     During  combustion  these 
two  agents  would  oppose  one  another  and  E  might  remain 
fairly  constant.     Experimental  studies  of  the  variations  of 
E  and  the  intrinsic  flame  speed  | F  -  u|    are  projected. 

Various  representations  for  p,  in  terms  of  time  integrall 
involTing  5,  Y  and  t  for  any  given  particle   (layer)  of 
the  loixture,  are  obtainable  from  (4)t 
If   Y   is  constant 

t 

(22)  P'^'^^t-O  "/  ^'^'^^  ^   '^^'^  ** 

0 

for  a  given  parti ole.  For  the  oase  of  non-burning 

(E  -  0)  at  time  t  (22)  yields  the  usual  adiabatio 

relation!  p  T  "^^  »  const.   If,  as  seems  to  be  a  good 

approximation,  C  »  R(C  +  G©),  where  R,C  and  G  are 

P 
constants,  then 


(22.)    ^  («)-=.  „p  (,.  »^)G  -/   ^« 

''-00 


with  pT  •  Rd  , 

where  p-.  lO-,  represent  the  values  of  p  and  9  outside 
the  reed  valves,  and  where  E  »  0  until  the  particle 
begins  to  burn. 

Relations  like  (22)  and  (22')  are  extensions  of 
pT    »  const,  to  cases  involving  combustion,  and  they 
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may  be  useful  if  approximate  values  for  the  integrands 
are  known.     For  example,  the  integrand  in  (22)  might 
be  regarded  as  constant  during  combustion  of  the 
particle,  and  zero  otherwise. 

(Iv)     Perturbation  aeries  for  the  equations.     When  the 
pulse- Jet   is   first   fired  the   combustion  disturbs  the 
initially  steadily  flowing   (or  resting)  mixture,   and  for 
a  certain  interval  the  changes  may  be  treated  as  fairly 
small  perturbations    (if  violent  detonations  are  avoided 
by  use  of  suitable  mixtures),      since  the  normal  combustion 
in  pulse-jets  appears  to  be  rather  mild  (S  small)   it  might 
be  expected  that  a  couple  of  terms  in  perturbation  series 
expansions^        of  the  dependent  variables  in  powers  of  a 
parameter   such  as  the   K  defined  below  would  provide   good 
approximations 

(23)  K  =  maximum  value  of  [(Y-  1)E]. 

It  is  also  convenient  to  introduce  I  =  (Y-  1)E/K  which 
therefore  has  the  maximum  value  unity.     The  perturbation 
series  in  powers  of  K  is  given  in  Appendix  II.      If 
sufficient   approximate  data  on  the  temperature  and  pressure 
distributions   in  a  pulse-jet   during  its   periodic   oscillations 
were  determined  experimentally  it  would  most  probably  be 
possible  to  deduce  a  mathematical  representation  which 
could  be  used  to  provide  much  better  first  terms  in 
perturbation  series  of  various  convenient  types.      Analytical 
rather  than  purely  numerical  approximations  are  preferred 
for  the  investigation  of  the  relative  importance  of  the 
various  physical  parameters  which  govern  the  phenomena  in 
pulse-jets. 

The  non-uniformity  in  the   cross-section  area  S(x)    of 
the  pulse- jet  leads  to  analytical  complications  which  are 
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mlnimlzed,   in  general,  by  treating  the  non-unifonnlty  as  a 
pertTorbat ion  of  the  simple  form   (a)   or   (b)   Indicated  in 
Pig.    6    (cf.    Fig.    3). 


CROSS- 
SECTIONS 


'So 

12 


(a)  Limiting  form  of  tube 
for  the  one -dimension- 
alized  pulse-Jet. 


(b)  Uniform  tube 
for  the 
pulse-jet. 


Pig,  6   Approximate  shape  of  pulse-jet  tube. 


In  using  the  form  indicated  in  Pig.  6(a)  as  the  limit 
of  a  smooth  one-diinensionalized  form,  continuity  of 
[s(x)u/t],  [u^  +  (2YpT)/(Y-  1)3  and  [pr*"]  across  the 
discontinuity  in  S(x),  is  implied  by  the  hydro-thermo- 
djniamical  equations  (5)-(7),  This  is  a  mathematical  rather 
than  a  strictly  physical  limit  (which  would  involve  three- 
dimensional  considerations  of  reflections  from  constrictions). 
The  mathematical  limit  probably  provides  a  good  approximation 
for  the  effects  due  to  the  fairly  gradual  constriction 
indicated  in  Pig.  3,  but  it  is  rather  difficult  to  obtain 
analytical  (rather  than  purely  numerical)  solutions  of  the 
hydro-thermodynamical  equations  by  use  of  this  limit.  The 
form  indicated  in  Pig.  6(b)  is  much  simpler  for  use  as  a 
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st  art  Ing  point  for  successive  approximations  to  include  the 
effects  due  to  the  mildly  non-uniform  cross-section  of 
standard  pulse- jets.   An  explicit  formulation  for  perturbation 
series  based  on  deviations  from  the  form  indicated  In 
Pig.  6(b)  Is  given  in  Appendix  II. 

Thus  a  double  perturbation  development  of  the  equations 
Is  proposed:   one  perturbation  being  due  to  the  effects  of 
combustion,  the  other  to  the  non-unlformlty  In  the  cross- 
section  of  the  tube.   In  the  simple  example  described  In 
the  following  paragraphs  only  the  first  of  these  perturbations 
Is  used. 

(v)  Explicit  treatment  of  an  Idealized  case.   An  example 
Is  Indicated  In  Fig.  7  belov/,  which,  in  view  of  observations 
related  to  Fig.  5,  probably  is  reasonably  Indicative  for 
the  main  characteristics  of  the  pulsating  V.l  motor. 


INTAKE   I  Pq^o^o'^o     i  ^   ^o'^l^o^l       EXHAUST 

x  =  0        x  =  a^  x  =  L 

o 

(point  of  re-lgnltion 
at  time  t  =  0) 


Pig.  7   Idealized  example  of  a  pulse -Jet 
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Speclficatlona ; 

TTnlform  tube  with  cross-section  area  S  and  length  L. 

Initial  conditions; 

At  time  t  =  0:  the  part  of  the  tube  between  x  =  0  and 
X  =  a   (the  point  of  re-ignition)  contains  a  fresh,  combustible, 
uniform  mixture  with  pressure  p  ,  specific  volume  T  , 
temperature  ©  ,  flow  velocity  u  ,  sound  speed  c  and  ratio 
of  specific  heats  Y  ;  the  rest  of  the  tube  (from  x  =  a 
to  X  =  L)  contains  a  hot,  completely  burnt  gas  with 
corresponding  physical  values  p  ,  t,,  6,,  u  ,  c,  and  Y,. 
These  initial  conditions  are  presumed  to  correspond 
roughly  to  the  conditions  in  an  actual  V.l  motor  at  the 
instant  of  re-ignition. 

Boundary  conditions; 

(cf.  pages  14  and  15):   At  the  exhaust  end,  x  =  L, 
p  =  p  .   At  the  intake  end,  x  =  0:  p  -  p.  =  -(u  -  u  )Z  if 
u  >  0  (open  valve  condition),  or  u  =  0  If  p  >  p  +  Zu 
(closed  valve  condition).   Here  Z  Is  a  constant  which  should 
be  determined  experimentally  as  indicated  in  Section  4, 
but  which,  on  the  basis  of  generalized  forms  of  the  simple 
Bernoulli  equation^   and  (19)  (page  15)  Is  of  the  ordei*  of 
u  /t  .  We  therefore  assume  z  =  "^^ "^q* 

Paeudo-Lagrange  coordinates  (m,t).   (See  Appendices  I  and  III) 

For  convenience  the  mass  flux  leading  to  the  above 

Initial  condition  may  be  regarded  as  corresponding  to  a 

constant  particle  velocity  u  throughout  the  tube  until 

t  =  0.  Taking  the  reference  time  t   (see  page  8)  to  be  the 

corresponding  time  of  arrival  (at  the  Intake)  of  the  particle 

which  reaches  x  =  L  at  t  =  0,  It  follows  that  t^  =  -L/u^ 

'  o     '  o 

and  that 


(L  +  u^t  -  x)s/'rT  If  X  >  a^  +  u^t 

0        '   1  CO 

(24)  m  = 

i^(L  -  ao)S/'r^  +  (a^  +  u^t  -  x)s/rQ  If  x  <  a^  +  u^t 
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untll  the  disturbance  due  to  the  burning  (which  starts  at 
X  =  a  at  t  =  O)  reaches  the  position  x   concerned.  The 
general  relation  connecting  in,  x   and  t  is  discussed  In 
Appendix  III, 

Combustion  and  flame  speed  conditions. 

The  combustion  term  I  defined  below  equation  (23) 
(page  19)  and  concerned  in  (7),  (8)  and  (11)  (page  11),  la 
here  assvuned  to  have  the  following  values: 

1  during  combustion  of  a  particle 

(25)  I  = 

0  before  and  after  combustion  of  a  particle. 

The  flame  speed  velocity  is  assumed  to  be  such  that  if 
m  =  ni(t)  is  the  pseudo-Lagrange  coordinate  of  the  flame 
front  in  the  fresh  mixture 

(26)  ^  m(t)  =  m  =  const. 

This  assumption  amounts  to  regarding  the  flame  spread  as 
analagous  to  the  successive  firing  of  fire-crackers  connected 
In  line  by  fuses  of  uniform  length.  Experimental  data  is 
required  to  determine  how  good  these  assumptions  may  be. 
After  t  =  0  until  the  flame  Is  extinguished  (cf.  Pig,  5, 
page  16)  the  burning  extends,  according  to  (26)  and  (24), 
from  the  particle  v;lth  m  =  (L  -  aQ)S/T]^  (  =  m^,  say)  to  the 
particle  with  m  =  m  +  m  t. 

Explicit  solutions  of  the  pseudo-Lagrangian  hydro-thermo- 
dynamlcal  equations  (11)  -  (13)  subject  to  the  above  conditions 
are  given  to  first  order  terms  In  Appendix  III,  for  the 
first  half  of  a  cycle.  The  calculations  are  not  carried 
further  since  physical  data  are  lacking. 


-24- 


Numerlcal  example.  Using  the  formulas  given  in  Appendix  III, 
and  the  follov/ing  typical  numerical  values  for  the  physical 
quantities  involved^  the  results  shown  in  the  graphs  In 
Figures  8  and  9  were  obtained: 

K  =  5.85'«10''"'''cin^/sec^,  Y  =  1.4,  c  =  3»10^cni.aec 

u_  =  .2c  ,  f  =  .3c   (based  on  an  analysis  of  Pig.  5, 
°      °        °  page  16) 

Z  =  V^o'  '^l/'^o  =  '^'  Y^  =  1.2 
c^  =  2.5c   (based  on  an  absolute  temperature  ©,  =  7©  ) 

a  =  45  cm  (1.5  ft),  L  =  300  cm  (lOft)  (cf.  Pig.  1,  page  l) 

Although  Fig.  8  corresponds  to  only  a  first  approximation 
for  propagation  phenomena  in  the  pulse- jet.  It  provides  a 
good  framework  for  the  analysis  of  actual  motions  in  the 
jet.   Thus  C,  and  Cg,  for  example,  represent  compression 
wave  fronts  which  are  reflected  from  the  open  exhaust  end 
as  rarefaction  wave  fronts  R,  and  Rg.   R,  is  reflected 
from  the  closed  end  as  a  rarefaction  wave-front,  follov/ing 
closely  behind  the  rarefaction  wave-front  R,  which  is  pro- 
duced when  the  flame  front  reaches  the  closed  intake  valves. 
These  rarefaction  wave-fronts  will  in  turn  be  reflected  from 
the  open  end  as  compression  waves.   Judging  from  Pig.  5, 
these  rarefaction  waves  appear  to  extinguish  the  flames  as 
they  pass  through  the  region  of  burning  (at  a  time  corresponding 
roughly  to  the  upper  limit  of  Pig.  8);  and  the  last 
mentioned  compression  waves  probably  stimulate  the  re-lgnltion 
automatically  at,  say,  t  =  20  ms.   Pig.  9  indicates  that  the 
valves  reopen  at  about  t  =  10.2  m.s. ,  allowing  a  fresh  gas 
mixture  to  enter  the  pulse- jet. 

The  right  side  of  the  flame  region  in  Pig.  8  Indicates 
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10  ma 


5  ms 


Pig.  8 


Diagram  for  the  first  approximation  to  the  wave  fronts 
(characteristics)  and  flame  region  relative  to  the 
pulse- jet  tube  of  Pig.  7.  The  first  disturbances,  due 
to  the  Ignition  at  x  =  1.5,  reach  the  Intake  valves 
at  time  t  =  1.9  ms,  and  the  exhaust  end  at  t  =  3.1  ms. 
The  Intake  valves  close  at  time  t  =  3.3  ms,  and  reopen 
at  time  t  =  10.2  ms.   The  flame  front  reaches  the 
Intake  at  t  =  7.2  ms;  the  shading  indicates  the  flame 
region.  The  Important  wave  fronts  are  Indicated  by 
solid  lines;  the  dotted  lines  represent  relatively 
weak  reflections  from  the  Interface  between  the 
previously  burnt  and  the  fresh  gases. 
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the  actual  path  of  the  interface  particles.  The  rapid 
acceleration  of  these  and  the  adjacent  exhaust  gases  is 
notable  (see  also  Pig.  9b)  and  probably  makes  u^  +  Kx^ 
'see  Appendix  III)  a  rather  poor  approximation  for  the 
gas  velocity  for  values  of  x  and  t  lying  above  the  double 
line  E  in  Pig.  8.   (Higher  approximations  can  be  calculated 
without  encountering  too  complicated  fomns,  but  have  not 
been  included  since  the  values  (and  the  nature)  of  some  of 
the  parameters  used  in  this  example  are  based  on  guesses 
rather  than  reliable  data.)  It  is,  for  example,  entirely 
unclear  Just  how  and  when  the  flame  is  extinguished. 
The  indicated  motion  of  the  flame  region  is  in  qualitative 
agreement  with  that  actually  observed  in  the  early  part 
(say  the  first  third)  of  the  cycle  shown  in  Pig.  5  for  the 
V. 1  motor.   For  a  comparison  between  theory  and  observation 
at  later  stages  of  the  cycle  it  would  appear  advisable  to 
modify  the  perturbation  scheme  at  some  point.   A  possible 
procedure  would  be  to  use  Rlemann*s^   exact  integral 
representations  of  solutions  of  (5)  and  (6)  subject  to  (3), 
starting  with  values  given  by  the  present  first  approximation 
on  the  above  mentioned  line  E  In  Pig.  8.   Graphical  and 
numerical  methods  are  also  available^   .  (See  App.  IV.) 

In  Pig.  9a  an  Interesting  feature  is  the  considerable 
difference  between  the  pressures  at  the  Intake  end  and  at 
the  interface  between  the  fresh  and  previously  burnt  mixture. 
This  difference  suggests  that  the  ass\imptlon  (used  In 
previous  theories^  '^) that  the  state  of  the  gas  is 
effectively  uniform  at  all  times  in  the  combustion  chamber 
Is  of  doubtful  validity. 

It  should  also  be  noted  that  the  pressure  curve  at  the 

Intake  end  for  the  first  half  of  the  cycle,  as  indicated  in 

(1,(2 


Pig.  9a,  is  in  satisfactory  agreement  with  rough  observations. 

In  Pig.  9a,  the  rapid  drop  in  p/p^  at  the  Intake  valves 
after  the  flame  front  arrives  is  notable.  Since  the  pressure 
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Fig.  9a,b. 


Graphs  for  pressures  and  velocities 
(first  approximation). 
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at  the  valves  Is  essentially  ^at  drives  the  pulse-Jet, 
this  phenomenon  suggests  that  the  thrust  might  be 
considerably  Increased  by  keeping  the  flame  away  from  the 
exhaust  end  for  a  longer  time  -  for  example  by  having  a 
lGa?ger  amount  of  combustible  gas  In  the  pulse-jet  at  the 
time  of  firing. 

The  rapid  acceleration  of  the  exhaust  gases,  mentioned 
above.  Is  Indicated  in  Fig.  9b. 

4.  Concerning  the  nature  of  desired  experimental  data  on 
pulse-Jets. 

In  order  to  make  a  reasonably  accurate  analysis  of 
pulse-Jets  in  terms  of  a  mathematical  framework  such  as 
that  indicated  in  the  preceding  paragraphs,  it  is  necessary 
to  secxire  more  fundamental  physical  data  than  is  available 
at  present.   It  is  hoped  that  investigations,  projected 
or  in  progress  at  N.R.L. ,  N.A.C.A. ,  Aerojet  Corporation  and 
N.Y.IJ.,  will  be  expanded  so  as  to  provide  the  desired  data. 
In  this  connection  the  following  comments  may  be  helpful: 
High  speed  movies  (2000  or  more  frames  per  second,  both 
direct  and  shadowgraphic  pictures):  Judging  from  the 
preliminary  movies  taken  through  windows  (or  glass  walls) 
in  pulse -Jets,  much  information  can  be  obtained  relatively 
easily  by  this  method.  Direct  high  speed  movies  can  show 
(quantitatively)  the  motions  of  the  flames  from  the  time 
of  re-ignltlon  until  the  time  of  extinction:  on  the  other 
hand,  if  a  sodium  compound  is  added  to  the  fuel  to  make 
the  hot  gas  glow,  a  measure  of  the  turbulent  motions  of  the 
hot  gases  may  be  obtained.   Shadowgraphic  high-speed  movies 
may  be  made  by  directing  sunlight  through  two  transparent 
sides  of  a  pulse-Jet  and  then  onto  an  otheirwise  shaded 
screen  which  can  be  photographed.   Such  shadowgraphic  movies 
can  show  the  turbulent  motions  of  the  gases  as  they  move 
from  the  intake  down  the  tube,  and  quantitative  information 
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may  thus  be  obtained  on  the  one-dlmenslonalized  velocities 
involved.  It  would  be  desirable  to  obtain  correlated 
direct  and  shadowgraphic  pictures  in  order  to  distinguish 
between  flame  speeds  and  particle  speeds,  and  it  would 
furthermore  be  desirable  to  use  a  mirror  system  to  permit 
the  viewing  (on  the  same  films)  of  the  motions  of  the  intake 
valves.  Ultimately  it  would  be  desirable  to  correlate  the 
above  observations  with  temperature  and  pressure  measure- 
ments.  (It  might  be  remarked  here  that  high  speed  movie 
studies  of  flame  propagation  in  tubes  under  controlled 
conditions  of  ignition  and  turbulence,  are  under  way  at 
N.Y.U. ) 

Intake  valve  conditions.  It  would  be  of  value  to  investigate 
experimentally  the  nature  of  the  boundary  condition  at  the 
reed  valves  and  to  see,  in  particular,  if  relation  (19) 
page  15  is  adequate.   It  should  be  relatively  simple  by 
means  of  Pltot  gauges,  to  determine  the  relation  among  the 
ram  pressure  and  the  velocity  and  pressure  just  inside  the 
intake  end  of  the  motor  under  steady  flow  conditions  in  a 
wind  tunnel  (without  combustion).  The  question  of  how 
pulsing  flow  might  modify  such  a  relation  because  of  dynamical 
factors  would  probably  be  fairly  hard  to  answer  experimentally. 
Possibly  the  dynamical  effects  might  be  evaluated  by  the 
use  of  shadowgram  movies  to  measure  effective  (variable) 
speeds  of  the  gas  (delineated  by  mass  movement  of  small 
vortices),  coupled  with  oscillograms  of  the  "static" 
pressures,  a  reciprocating  pump  being  employed  to  produce 
an  "artificially"  pulsed  (cool)  exhaust. 
Pressure  and  temperature  oscillograms.   Standard  pick-ups 
for  measuring  variable  ("static")  pressures  and  variable 
temperatures  in  gases  are  apparently  subject  to  large 
perturbations  under  the  conditions  of  large  and  rapid 
mechanical  and  thermal  oscillations  encountered  in  pulse- 
jets.   An  important  problem  is  therefore  to  design,  constmact 
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and  calibrate  suitable  measuring  equipment,   some  progress 
along  these  lines  has  been  made,  in  particular  by  Aerojet 
and  by  the  AMG-NYU  in  collaboration  with  Ensign  Sdelman 
of  BuAer.  Various  methods  for  avoiding  the  effects  of 
'^su'iation  of  the  elastic  "constants"  with  temperature  in 
pressure  diaphragms,  have  been  incorporated  in  the  designs. 
Photo-electric  scanning  is  now  being  investigated  in 
connection  with  spectroscopic  methods  for  recording  variable 
temperatures  in  hot  gases.   After  such  measuring  devices 
have  been  precisely  calibrated  under  conditions  simulating 
those  in  pulse-jets  It  would  "be  desirable  to  secure 
oscillograms  of  the  pressures  and  temperatures  simultaneously 
at  several  representative  positions  in  pulse-Jet  motors. 
This  information  would  permit,  in  particular,  determination 
of  the  energy  release  rate  fe  (cf.  equations  (4),  (7),  (8) 
and  (11))  and  of  variables,  other  than  velocities,  involved 
in  initial  and  end  conditions.   (More  convenient  than 
oscillograms  would  be  movies  of  strings  of  glow  lamps 
suspended  above  the  motor  and  In  view  of  the  high  speed 
cameras  being  used  as  described  In  a  preceding  paragraph. 
Circuits  can  be  arranged  to  make  the  number  of  lamps  glowing 
in  a  string  at  each  Instant  proportional  to  the  pressure 
or  temperature  signal  concerned. ) 

General  data  required  for  comparison  of  theory  and 
actual  performance; 

(1)  Dimensions  of  components  of  motor.   (11)  Constants 
Involved  in  Intake  (reed  valve)  boundary  condition,   (ill) 
Rate  of  flow  of  fuel  (preferably  oscillogram  record). 
(iv)  Combustion  and  flame  constants  for  the  fuel  under  the 
turbulence  conditions  involved  (determined  by  pressure, 
temperature,  particle  velocity  (high  speed  movies)  and 
flame  pattern  observations  in  pulse- jets  in  general). 
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APPENDIX  I 


Derivation  of  the  equations  In  (in,t), 


(1)  At  time  t,  the  mass  contained  between  the 


z' ■ 

I 'STREAM 
I  TUBE'  OP 
I  GASES 


m 


m+4!n 


X  x+/ix 


layers  at  x  and  at  x  +    ax  is  -   Am.     Hence 


-   Am  =  S(x)»  Ax/r  • 


Thun,  in  the  limit,  we  have  the  continuity  equation 


(12) 


X  =  -r/s(x). 
m 


(ii)  The  equation  of  motion.  Force 
for  a  small  layer  of  gas  is 


=  mass  X  acceleration. 


S(x)[p(m)  -  p(m  +Am)]  =  -(Am)x.  .. 


Hence,  in  the  limit. 


(13) 


S(^)Pm  =  ^tt' 
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(111)  Since  m   Is  constant  for  a  given  particle  (layer), 
the  internal  energy  equation 

(11)  (Y-  1)E  =  Ypr^  +  -rp^ 

remains  formally  the  same  as  in  Lagrangian  coordinates. 
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APPEMDIX  II 
Perturbation  Series 


,S(x)-S  , 
1.   Perturbation  series  in  e  =  maxl  ■■■-■r  -r  -■  • 

The  constant  S  is  some  conveniently  chosen  mean  value  of 
S(x).  We  suppose  that  the  cross-section  of  the  pulse-Jet 
tube  does  not  vary  too  much,  and  hence  that  e  is  a  small 
number.  For  the  German  VI  motor,  e  ~.3,   A  power  series 
development  of  x,  p,  and  t  in  terms  of  e  is  then  assumed: 


X  =  X  +  eX  +  enc   +  ... 

(26)  p  =  P  +  ep*  +  e^p**  +  ... 

™*    2™*-*  . 

T  =  T  +  eT  +  e^   +  ... 

where  all  quantities  except  the  small  constant  e  are 

functio 

30  that 


functions  of  m  and  t.  Furthermore  we  may  define  S  (x) 


S(x)  =  S^  +  es*(x). 


Then 


(27) 


S(x)  =  S^  +  eS*(X  +  ex*  +  ...  ) 


o 


=  Sq  +  es*(x)  +  e^X*i4(x)  + 


We  now  substitute  (26)  and  (27)  into  (11),  (12)  and  (13) 
and  equate  coefficients  of  like  powers  of  e.  This  furnishes 
the  following  sets  of  equations. 
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^o^ 


•^  =  -  ^o^ 
tt  o  m 


(Y-  1)E  =     YPT^  +  TP^ 


(terms  not 
involving  e) 


En 


At  =  s*(x)p„  +  S^P^ 
^   Y(PT*   +  P%)   +   (TP*   +  T*P^) 


=  0 


(coefficients 
of  e) 


T        =  -[SqX 


^2    ^ 


_•» 


-*_* 


-** 


_*«• 


(coefficients 
of  e^) 


Y(PT.      +  P  T^  +  P     T^)   +   (TP^     +  T  P^   +  T     P^)   =  0   ,   etc. 


We  note  that  E     is  a  non-linear  system;     but  once  E 
Is   solved  so  that  X,  T,  and  P  are  known  functions  of  m  and  t, 
the  successive   systems  E, ,  Ep,  etc.   are  each  linear  and  may 

it  it  it  itit  itit  itit 

he  solved  in  turn  for  X  ,  P  ,  T*;  iC    ,   P  ,  T   etc.  The 
system  E  corresponds  to  a  pulse-jet  of  uniform  cross-section 
S  .  Since  it  is  non-linear  we  seek  a  second  perturbation 
series  for  its  solution. 

2,   perturbation  series  for  E  in  K  =  max(Y-  l)E. 

We  suppose  that  a  power  series  development  of  X,  P  etc. 
in  powers  of  K  is  possible  so  that 


=  %+   K%*  +  K^%**  +  .. 


(28) 


X 

P  =  (?+  K  (?*  +  ..  . 

T 

Y  =  Jb  +  KiJ"  +  ... 


=  7  +  K7   +  ... 
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We  define  I  by  the  equation 


so  that 


I  =  (  Y  -  1)E/K  , 


(29) 


1   =  d    +   KcX*(P,T,t) 

=  c8-  +  K^*((?,7,t)  +  I?ie*d*   +  7*d*)   +  ... 


and  c?.  =  1  or  0  according  as  the  particle  is  burning  or  not. 
Proceeding  as  before,  we  substitute  (28)  and  (29)  into  E  ,  p 
and  equate  coefficients  of  like  powers  of  K.  We  obtain 


K 


r 


o  n 


V 


^t  =  s,P^ 


^(?^   =  0 


(terms  not  involving  K) 


Ei 


tt 


% 


o  m 
o  n 


(coefficients  of  K) 


S.    = 


\ 


.^*((?Jt^  ■"  ^^t  •"  '^*^t 


etc. 


As  before,  only  E*  is  non-linear,  the  succeeding 


systems  being  all  linear.  The  system  E'  corresponds  to  a 
pulse-Jet  of  uniform  cross-section  with  no  burning  —  i.e. 
is  merely  a  set  of  relations  which  the  initial  state  of  the 
system,  given  by  (?,  CT  and  %,   must  satisfy  at  the  time 
t  =  0  of  firing.  We  are  hence  not  required  to  solve  E* ,  but 
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merely  to  check  that  the  initial  conditions  (i.e.  the 
functions  <P  ,  J  and  %  )  which  we  postulate,  are  in 
agreement  with  these  relations. 
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APPENDIX  III 
Explicit  Solution  of  an  Idealized  Case 

1«   General  Solution. 

We  consider  the  idealized  pulse-Jet,  described  on 
pages  21  and  22,  with  uniform  cross-section  S  and  length  L. 


x=0         a^  L 

(point  of  reignitlon 
at  time  t  =  O) 


Pig.  10   Idealized  pulse-Jet, 


At  the  time  t  =  0  of  firing,  it  is  assvuned  that  the  gas  In 

the  tube  is  In  motion  with  constant  velocity  u  and  has 

o 

constant  pressure  p  ;  and  that  a  fresh  combustible  mixture 
with  constant  specific  volume  t  and  adiabatic  exponent  V 
extends  from  x  =  0  to  x  =  a  ,  the  point  of  firing,  while 
completely  burnt  gas  with  constant  specific  volume   T,  and 
adiabatic  exponent   Vt  extends  from  x  =  a  to  x  =  L«   Since 
the  cross-section  is  uniform,  S(x)  =  S  =  S,  and  e  =  0. 
Hence,  referring  to  (26),  x  =  X,  p  =  P,  and  t  =  T.  According- 
ly, only  the  perturbation  series  In  powers  of  K  is  to  be 
considered.   (In  considering  a  pulse-Jet  of  non-unifonn 
cross-section,  the  present  solution  would  of  course  be  the 
starting  point  for  the  perturbation  series  in  e.) 

Prom  the  prescribed  initial  conditions  we  have 
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7  = 


0  <  X  <  a^ 

a^  <  X  <  L 


(?=   p. 


>^  = 


V. 


0  <  X  <  a 
a^  <  X  <  L  . 


We  also  have,  as  is  clear  from  inspection  of  Pig.  11, 


m 


o 

m 


H 


Pig,  11   Relation  between  m  and  x  at  time  t. 


''  (L  +  u^t  -  x)s/t^ 


m  =   < 


a^  <  X  <  L 
o  —   — 


""O  •*•  ^%  -^  %*  -  '^^S/T^     ^1^1%    f 


where  m     is  the  mass   of  burnt  gas. 


m 


=   (L  -  &^)s/r^ 


Hence 


%   = 


0  <  m  <  m^ 
—      —     o 


L  +  u^t   -  mT,/s 
o  1' 

y%  +  u^t   -   (m  -  m^)  T^S  m  >  m^ 
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We  observe  that  V ,  (?,  M    and  %  so  prescribed  satisfy 

the  system  E' . 
o 

bince  we  are  taking  m  and  t  as  independent  variables, 
we  must  refer  the  boundary  condition  at  x  =  0  and  x  =  L 

to  the  corresponding  points  in  the  (in,t)  plane.   Now 
x=  %  +  KX      +...  ,  and  it  follows  that  the  value  of  m 
for  the  particle  which  reaches  x  =  L  at  time  t  will  be  given 

m   =  u^St/ r^  +  [terms  involving  K] . 


For  the  first  approximation,  then,  we  replace  the  actual 
boundary  x  =  L  by  the  approximation 


m  =  U^St/ T  T 
o  '   1 


Similarly,  for  the  boundary  x  =  0  (the  intake  valve  end) 
we  take 

m*  +  u  St/  T-      until  the  valves  close 
o    o  '^   1 

1  mg  after  the  valves  close 

where  m'  is  the  total  mass  of  gas  in  the  tube  at  time  t  =  0, 


m 


;=  [(L  -  a^)/T^  +  V-^o^S  , 


and 


"^h^K""   ^0^4/  '^o  ' 
tX   being  the  time  of  closure  of  the  valves. 
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We  assume  that  the  flame-front  is  propagated  from 
particle  to  particle  at  a  constant  rate  f  —  Independent 
of  the  motion  or  density  of  the  gas.   If  m-(t)  is  the  value 
of  m  for  the  parricle  which  begins  burning  at  time  t,  we  have 


m-.(t)  =  m^  +  fSt/  T^  =  m^  + 
r      o     '00 


mt 


with  m  =  fS/  T^. 

In  the  (m,t)  plane,  then,  the  region  in  which  we  wish 
to  solve  the  system  Ei  is  as  indicated  in  Pig.  12  (below). 

The  system  Ei  can  now  be  written  in  specific  form, 
including  the  initial  conditions: 


/' 


^tt  =  sc?: 

'^*  =  -  s%* 

s»    « 

m 

1 

roPo't  *%Pt=^ 

m^  <  m 

0 

.  Vo^t  *'fiPt  =  ° 

m  <  m 

where,  at  any  time  t. 


c^  = 


m     <  m  <  m- 
o  f 


in  >  m^  . 


The  first  equation  in  E*   can  be  satisfied  identically 
by  introducing  a  function  V(m,t),   analogous  to  a  stream- 
function,   such  that 


%.    =  const   •    S  V„ 
b  in 

(P     =  const   •   V. 


t 


-41- 


Oil 

4>| 


II 

El 


Burnt  Gas 


o 
o 


Flame  Front 
Reaches  Valves 


Burning  Gas 


'i' 


Valves  close  t 


Fresh  Combustible  Gas/^ 


m. 


Fig,  12   Region  In  (m,t)  plane  corresponding  to  pulse-Jet. 
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It  is  convenient  to  choose  the  constant  so  that 


o^     ' 


*      1 
o 


where  a  is  given  "by 


«^  =■  Vo^^/  ■^o  • 


In  terms  of  V,  we  have 

,2 


3t  =  -  S'^.t  =  2.  (,S.i2,  \  ' 


the  system  Ei  then  reduces  to  the  pair  of  equations 
2 


-YiP-S 


2 


S^  ^tnm  ^  o».  r}   '5,  ^^t  =  0    ^  <  ni^  • 


2^^{&''-m^)     "^       2T^(a''-ni^) 


These  may  he  conveniently  written  as 

(31)  \^   -  a^^  =  2{e^'m^)  ^  m  >  m^ 


(32)  V. .  -  b^V^„  =0  m  <  m, 

zz  mm  ( 
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where 


b^  =  Y,p,s2/  T. 


'1^0' 

Since  the  sound  speed  c  is  given  by 

c^  =  Xp/^      =  YpT  , 

we  have 

a  =  c^S/  r^   ,  \>  =   c^^S/  ^^  , 

c  and  c,  being  the  sound  speeds  in  the  fresh  and  burnt 
mixtures  respectively. 

Equations  (31)  and  (32)  are  hyperbolic  partial 
differential  equations  whose  general  solution  may  be  written 
in  the  form 


(33)   V  =- 


-  •   2 

f(m  +  at)  +  g(m  -  at)  -  (m  -  m  -  mt)  (9^  m  >  m 

^{m   +  bt)  +  G(m  -  bt)  m  <  m 


•   2 

The  function  -  (m  -  m  -  mt)  -J^  (  ^  defined  on  page  40 )  la 

a  particular  integral  of  (31),  which  vanishes  together  with 

its  first  partial  derivatives  along  the  flame  front 

n  -  m  -  mt  =  0. 
o 

The  first  approximation  for  the  pressure  and  velocity 
may  now  be  written  in  terms  of  the  function  V, 


'=*'  ^  =  "o  ^  ^:^p^  '* 


(35)  U  =  U^  +  ^n   H  V   , 
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It  remains  only  to  determine  the  function  V  (i.e.    f,   g,  P 
and  G)   in  such  a  way  as  to  satisfy  the  boundary  conditions. 

Boundary  Conditions. 

(l)      Open  end.      At  the  open  end,  we   should  have 

P  =  Po>   i«©- 

(36)  ^t  ~  ^       at  m  =  u^St/  T^  • 

(ii)     Intake  end.      At  the  intake  valves,  we   should  have 
p  -  p     =  -2(u  -  u   )  =  -Uq(u  -  u^)/  r^,  until  the  valves 
close   (at  t  =  tkt  when  u  =  O);     and  thereafter  u  =  0.     Hence 

(37)  V^  =  -  Su^vy  r^  at  m  =  m^  +  u^St/  r^   ,   t  <  t^ 

(38)  Vj^  =  -  2  TQUQ(a^  -  m^)/KS       at  m  =  m^  ,  t  >  t^ 

We  further  require  that  V,  V  and  V.  be  continuous 

throughout  —  and  in  particular  at  the  interface,  m  =  m  . 

At  t  =  0;  V,  V  and  V.  are  zero, 
m     z 

The  characteristics  for  equations  (31)  and  (32)  are 
straight  lines  of  the  form 

m  +  at  =  const.     m  >  m^ 

-  o 

m  +  bt  =  const.     m  <  m^  , 

—  o 

In  Pig.  13  below,  we  have  redrawn  Pig.  12,  but  included  all 
relevant  characteristics. 

Referring  to  Pig.  13  and  (33),  it  may  be  shown  that  V 
will  have  the  form 
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Pig.    13       Region  in   (in,t;   plane  corresponding  to  pulse-jet 
including  relevant  characteristics,    Cf,   Pig.    12. 
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(39)   V  =  ^ 


\ 


-  c?..«(in  -  in  -  mt)   +  -C(m  -  ra  -  ab)   + 

^(m  -  m  +  at  -  at^)   +  ...  m  >  m 

B*(m  -  m^  +  bt^   +  Bgfn  -  m^  -  bt  +  bt^)^  +  ...   m  <  m. 


where  the  A  's  and  b''  's  are  zero  for  values  of  m  and  t  (in  Pl^^.  13) 

below  the  corresponding  characteristic  line,  and  have  constant 

values  A^jApjB-jBg  etc.  above  the  corresnonding  line.   This 

form  is  suggested  by  the  follov;ing  facts: 

(i)  that  V,  V  and  V.  must  be  continuous  across  all 
'  m      t 

characteristics,   (ii)  tho  particular  solution  — 
— (m  -  a  -  mt)  'cJ?  is  quadratic  in  m  and  t,  and  all  terms 
must  match  along  the  boundaries.   The  solution  is  of  course 
unique,  so  that  all  that  need  be  done  to  justify  (39)  is  to 
verify  that  the  constant- A' s  and  B's  can  be  determined  in 
such  a  way  as  to  satisfy  the  boundary  conditions.   The 
quadratic  form  may  also  be  deduced  directly  by  use  of  (33), 
but  this  is  rather  tedious;  we  will  assume  the  form 
Indicated  in  (39)  and  proceed  to  the  determination  of  the 
constants. 

Determination  of  A' s  and  B's. 

(i)   Interface.   Along  m  =  m  ,  V  and  its  derivatives 
are  required  to  be  continuous.  V/e  have 

p  •   2   n 

A,*(m  -  m  -  at)   —  (m  -  m  -  rat)  •«  m  >  m 


V  = 


Bv(m  -  iti  +  bt)  m  <  m^ 

1     o  o 


Hence  setting  m  =  m  and  remembering  that  cV  =  1  along 
m  =  ni_,  we  have 
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A  o2.  2        'S.  2  _  „  .2.  2 


A^at   -  mt  =  B,bt   • 
We  may  solve  for  A^  and  B- ,    obtaining 

A^  =  m(b  +  m)/a(a  +  b) 
Bj^  =  m(a  -  m)/b(a  +  b) 


Similarly, 


^4  =  ^  ^  =  ^1^ 

2a2 
^3  =  b(a+b)   ^  =  ^2^ 


A5  =  M^Ag 


B5  =  M2A3 
2 


^^  -  a(a+b)  ^2 
n.0       F+a  *6  • 


(li)     Exhaust  end.      At  the  exhaust  end,   just  after 
the  first  wave  arrives  at  t  =  fc?   ,  we  have 

o  2 

V  =  B,«(ni  -  Kq  +  fet)      +  B2«(in  -  m^  -  bt   +  btg) 

=  B^.[(ni  -  mj  +  b(t   -  tj)]2  +  g^.f^  .  ^„   .  ^,(^   .  ^„j32   ^ 
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aubstltutlng  In  (36),  we  obtain 

V^  =  2bCB^«(b  +  u^S/  -r^)  +  Bg.Cb  -  u^S/  T^)](t  -  tj)  =  0 


Hence 


b+U_S/T, 

B  =  2 i  Bt  =  -  M,B.  . 

2   b-u^S/r,  ^     3^ 


Similarly, 


B^ M3B3  ;  Bg MgB^  . 


(ill)  Intake  valves.  In  the  neighborhood  of  the  Intake 
valves,  just  after  the  first  disturbance  arrives  at  the 
valves  but  before  they  close,  V  is  given  by 

2  \2 

V  =  A^« (m  -  m  -  at)  +  Ag* (m  +  at  -  m^  -  at^) 

=  Aj^.Cdn  -  mj  -  a(t  -  tj^)]^  +  Ag'Cm  -  m^  +  a(t  -  t^)]^ 
Substituting  in  (37),  we  obtain 

2a[-  A^*(UqS/  t^  -  a)  +  Ag.  (u^S/  t^  +  a)](t  -  t{) 
-2Su. 


o 


'^o 


[A3^.(u^s/Tq  -  a)   +  A2.(u^S/t^  +  a)](t  -  t{). 


Hence, 


«  a+u^S/r^         ■"• 

o  '    o 
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After  the  valves  close,  at  (mAf  tA)  we  have 

V  =  A^.(m  -  n^  -  at)^  +  A^.dn  -  m^  +  a(t  -  t^))^ 
+  A^'Edn  -  n^  +  a(t  -  t^)]^  . 

Substituting  In  (38),  we  have 

V^  =  const  +  2a[-  A^  +  Ag  +  Ag](t  -  t^)   =  const. 

Hence 

Ag  =  A^   -   Ag   . 

Similarly,  we  obtain 

This  completes  the  determination  of  the  A's  and  B's. 
Determination  of  the  t*s. 


The  time  tg  corresponding  to  the  closing  of  the  valves 
jtermined  b^ 
have  from  (38), 


is  determined  by  the  equation  u  =  0  at  m  =  m*  +  u  St/  t  .  We 


V  ■  -  2  T  u^(a^  -  m^)/KS  , 
m       0  0        '/         w 


Hence 

T„u^ ( a^-m^ )  ( a+b )  ( a+u. S/t  ) 
2  "  ^      2KSm(b+m)(a-u^S/T^) 
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Thls  may  be  written  in  terms  of  the  quantities  f,  c  ,  and 
c,  as 

time  of  closure  of  valves 

^,  ^  ^,  ^  %<%-^  ^(^^l^^l'^o^^^-*-^^  ^ 
2    1     2Kf(c^Vfr^)(c^-u^) 

Ml  other  points  (m,t)  indicated  in  Fig.  13  are  determined 
oy  elementary  analytic  geometry,  as  the  intersections  of 
given  lines.  We  have: 

(i)  Intake  valves: 
first  disturbances  arrive 


o-s 


^1   c  -u.  » 


flame-front  arrives 


a.+u  tJL 
T 


^f   =     o     o  2 


rarefaction  from  tail-pipe 

a  +u  tA 

t»  =  t„  +  Q  ^Q  2 
6    3      C^ 
O 

miscellaneous  reflections 


a  +u  t* 
o 

4.t  -  4-    .    0    0  2 

*5    *2  ^  ^n 


(ii)     Flame  front: 
miscellaneous  reflections  cross  flame-front 


at,  at„  at,  at. 

•         »   ^2  ~  •         »   ^3  ~ 

m+a  m+a  a-m 


4. til  _  'J-    ,   i.tti  _  2    ,   i.m  —  1    .   i-ttt  _  4 

1  "  i.      '      2  ~  . _    '   ^3  ~  ^  1   '   ^4  ~  a+m 
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(lii)  Interface: 
mlscelleneo-gg  reflections  cross  interface 


a  +u.tA 
3  -  ""1  '  ^4  "  ^1  "  "  "c" 


t,  =  2t"  ;  t-  =  t,  +  2  °  °  ^ 


o 

(iv)  Exhaust  end: 
first  disturbance  arrives 

L-a 


1    C^+U^ 


subsequent  reflections 


C,t,  C,tr> 

^2  ~  ^1  ^  cThT:  '  ^3  -  ^3  ^  c7+ir  • 

10  1   o 


We  have  thus  completely  determined  the  function  V 
for  the  portion  of  the  cycle  corresponding  to  Fig.  13, 
Prom  this  function  V,  may  be  deduced  the  pressure,  position, 
etc.  of  any  particle  at  any  time.   In  particular,  we  are 
interested  In  the  following  functions  of  the  time 

(l)  the  velocity  of  gases  at  the  exhaust  end 

(li)  the  pressure  at  the  Interface 

(ill)  the  pressure  at  the  Intake  valves 

(Iv)  the  position  of  the  interface 

(v)  the  position  of  the  flame-front. 
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The  first  three  of  these  quantities  follow  immediately 
from  (33)  and  (34),  the  relations  between  the  A*s  and  B's 
and  the  formulation  of  V.  We  have  (with  each  term  being 
added  to  the  terms  above  it  after  the  time  indicated  in 
the  inequality) 


/ 


t  >  0 


+  B3^(t 


(i)  u 


2KS(b+u  S/t, ) 
=  u.  +  rr-^ — =-<   + 


Va^Ii^) 


B3(t 


+  B5(t 


-  tj)  ,  t  >  tj 

-  tS)  ,  t  >  t" 


-  t")  ,  t  >  t" 


\ 


etc. 


(ii)  J£-  =  1  +   o 


c-^(a^-i^) 


+  (A^  -  i^/a^)t 


4 


\ 


t  >  0 


+  (Ag  +  A^)(t  -  t-^)  t  >  t^ 

+  (A3  +  A5)(t  -  tg)  t  >  tg 

+     A^it      -   tg)  t   >  tg 

+  (Ag  +  A^Q)(t  -  t^)  t  >  t^ 
etc. 
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^         0  t  >  0 


(iii)  I- 


=  1  + 


Y.Ka' 


'o 

"T' u  '  n 


+  A3(t  -  t^) 


+  2A^(t  -  t^) 


(Aq  -  m^/a^lCt  -  t|) 


+  SAgCt  -  t^) 
+  2A^(t  -  t») 
etc. 


t  >  t^ 


t  >  tj 


t  >  t! 
4 


t  >  t» 


t  >  t| 


The  position  x(ni,t)  at  time  t  of  a  particle,  which 
corresponds  to  the  mass -coordinate  m,  is  given  by 

nt 
x(m,t)  =  x(m,0)  +  /   u(ffi,t)dt  . 

Using  (35),  this  reduces  to 


(39) 


x(m,t)   =  x(m,0)   +  U-t   +  ^^n — rr- 

°  2t  (a^-m^) 


f  \ 


dt 


In  particular,  for  the  interface,  m  =  m  and  x(m,0)  =  a  , 
so  that  (39)  becomes 
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KSa 


(IV)  x(..t)  =  a„«„t  .  ^^  ^^S_.ii, 


/ 


< 


v.* 


+   (Ag-J4o)(t-t4)^ 
etc. 


t  >  0 

t  >  t, 

t  >  t, 

t  >  tj 

t  >  t. 


The  mass-coordinate  for  the  particle  which  is   at  the  flame- 
front   at  time  t   is  given  by 


m  =  m^  +  mt   . 
o 


Substituting  this  value  for  m  in  (39)  we  obtain,   for  the 
position  of  the  flame  front,  the  expression 


(v) 

x(m,t) 


=  a  +(u  -f)t  + 


KS 


2Ta(a  -m   ) 


^  -  A^(m-a)^t^ 

+  A2(m+a)2(t-tJ'  f 
+  A3(m+a)^(t-tg»  f 


■°   ^        "   '        -  A4(ft-a)^(t-t^«  f 
^+  Ag(i+a)2(t-t^'  f 


t  >   0 
t  >  tj* 
t>t- 
t>  t- 
t  >  tj' 


2,   Wumerical  Example. 

Using  the  values  for  the  parameters  listed  on  page  24, 
we  obtain  the  following: 
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(i)  Velocity  of  gases  at  tail  pipe;   in  cm/ms  : 


u  =  6  +  17 


< 


1 

0 

t   >   0 

+    .43(t    -   3.1) 

t   >   3.1 

-    .27(t   -   6.6) 

t  >  6.6 

^   +   .86(t   -  8.1) 

t  >  8.1 

(ms.) 


(ii)   Pressure  ratio  p/p^  at  interface: 


r 


^  =  1  + 


.056t 

-  .035(t  -  3.7) 
-i       +  .llO(t  -  5.5) 

-  .096(t  -  6.2) 
+  .016(t  -  8.1) 


t  >  0 
t  >  3.7 
t  >  5,5 
t  >  6.2 
t  >  8.1 


(ms.) 


(ill)  Pressure  ratio  p/p-  at  intake  valves: 


^=  1  + 
P« 


r  0 

+  .039(t  -  1.9) 

+  .211(t  -  3.3) 

<      +  .060(t  -  5.9) 

-  .390(t  -  7.2) 

-  .200(t  -  7.6) 
^  -  .192(t  -  8.4) 


t  >  0 
t  >  1.9 
t  >  3.3 
t  >  5.9 
t  >  7.2 
t  >  7.6 
t  >  8.4 


(ms.) 
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(iv)  Position  of  the  interface  (cm) 

/^   .154t^ 

-  .095(t  -  5.7)' 
X  =  45  +  6t  +  .357  <(  +  .321(t  -  5,5)' 

+  .096(t  -  6.2)' 

^  +  .044(t  -  8.1)' 

(v)  Position  of  the  flame  front  (cm) 

/  -  .073t^ 


X  =  45  -  3t  +  .357  \ 


-  .110(t  -  2.9)' 
+  .357(t  -  4.2)' 

-  ,0l5(t  -  5.3)' 
^  +  .051(t  -  6.2)' 


t  >  0 
t  >  3.7 
t  >  5.5 
t  >  6.2 
t  >  8.1 


(ms,) 


t  >  0 
t  >  2.9 
t  >  4.2 
t  >  5.3 
t  >  6.2 


(ma.) 


These  are  the  formulas  used  in  drawing  the  graphs  in  Pig.  9 
and  the  burning  region  in  Pig.  8. 
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APPENDIX  IV 

On  Analytical  and  on  Numerical  Methods  for  Direct  Treatment 
of  the  One-dimenslonallzed  Thermo-hydrodynamlcal 
Non-Linear  ISquations 


As  was  mentioned  in  previous  sections,  analytical 
rather  than  pnrely  numerical  solutions  of  the  equations  are 
desirable  in  order  to  evaluate  simply  the  relative  importance 
of  the  various  physical  parameters  involved  in  pulse-Jets. 
Unfortunately,  explicit  analytical  solutions  of  the  non- 
linear equations  are  obtainable  at  present  only  in  rather 
special  cases.  For  any  explicit  numerical  case  it  Is 
possible  to  apply  numerical  difference  equation  (or 
graphical)  methods,  even  if  empirical  curves  are  Involved. 
Whatever  the  method  it  apparently  must  involve  considerations 
of  the  propagation  of  discontinuity  fronts  or  'characteristics' 
such  as  are  indicated  below. 

For  the  one-dimensionalized  pulse-Jet  the  thermo- 
hydrodynamlcal  equations  (for  a  not  necessarily  perfect  gas) 
take  the  following  forms  in  terms  of  Bulerlan  coordinates 
(x,t)  and  dependent  variables  p, T , u (cf .  p.  10,11): 

1^  ©  =  H  (p,r)  (equation  of  state)  (©  =  temperature). 

2.  E  -pDr  =  DU  (equation  of  internal  energy  IT) 

^^  ~  ^  "*"  ^"^  ^^   Eulerians,  or  D  =  sr 
in  pseudo-Lagrangians  (m,t)).   If  the 
internal  energy  U  (per  unit  mass)  is 
expressed  in  terras  of  p  and  t  then 
DU  =  U  Dp  +  U«-Dr  . 

3.  T.  +  wr    -  Tu_.  -ruti*  (x)/s(x)  =  O  (equation  of  continuity). 

4.  u.  +  uu_  +  Tp_.  =  0  (equation  of  motion). 
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5.  P  =  u  +  f  (equation  of  flame  speed  P  relative  to  the 

jet  tube,  in  terms  of  the  'intrinsic  flame 

speed*  f  relative  to  the  gas. 
In  the  above  equations  E  and  f  depend  on  local  turbulence 
and  other  parameters  which  require  experimental  study. 
Equation  i.  depends  on  the  composition  of  the  gas.   If  there 
are  no  shocks  p  and  u  are  continuous,  and  if  there  are  also 
no  temperature  discontinuities  (such  as  occur  at  the  inter- 
face between  fresh  and  previously  burned  mixtures)  then  f 
is  continuous.  Derivatives  of  such  continuous  functions  may 
Involve  discontinuities  across  certain  curves  (* characteristics M 
whose  differential  equations  In  terms  of  dx  and  dt  may  be 
found  as  illustrated  below: 
In  the  customary  notation  let 

6.  H]  =  jump  in  the  value  of  a  function  H  across  a  discontinuity. 

Then  if  p,u,T  are  continuous  across  a  curve  whose 
direction  in  the  (x,t)  plane  is  determined  by  the 
ratio  of  dt  to  dx  : 

7.  p  ]dx  +  p^]dt  =  dp]  =  0  ,  u^]dx  +  u^]dt  =  du3  =  0, 

T  ]dx  +  T.  3dt  =  dr]  =0  (assuming  that  P^'^t'^x'^t' 

f  f   and  T.  are  continuous 

X  u 

on  each  side  of  the  discontinuity 
curve. ) 
Furthermore  by  2,-4. 

8.  Ht^]   +     T^lvi]    (p  +  Vr)   +    \Pt^    +  Px^^l^p  -  ^^ 
+     r]u)     -    u^]T    =  s»(x)] 'ru/S(x) 


'^M 


p^]T     =   0 


The  relations  8.  may  be  used  to  eliminate  any  three  of  the 

six  jumps  Pt^»Px^»'^t-'''^x-'»"t-'»^x-'  ^"  '^'   ^°^  ^^^^  *°  provide 
three  linear  equations  for  the  remaining  three  jumps.  The 

three  equations  are  homogeneous  if  E]  and  S*(x)]  are  both 
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zero,  which  is  true  if  cbc  ^  Fdt  and  dx  ^   0  (i.e.  if  the 

characteristic  is  not  a  flame  boundary  or  is  not  a  fixed 

X  where  the  cross-soction  has  a  discontinuity  in  its  slope). 

If  the  equations  are  homogeneous  the  determinant  of  the 

system  must  vanish  for  at  least  one  of  the  jumps  to  he  non- 
zero.  This  determinant  relation  yields  the  following 
characteristics  (in  differential  form). 

9.  dx  =  udt   (i.e.  dm  =  0  in  pseudo-Lagrangian  coordinates), 

10.  dx  =  (u  +  c)dt  with  c  =  r/  (p+U^)7u   the  local  sound 

speed  relative  to  the  mediam. 

For  9.  the  relations  7.  and  8., with  S]  =  S'(x)]  =0, 
make  all  the  jumps  zei*o  except  possibly  t  ]  and  T.  ]  (but 
T^]  +  T^lu  =  0). 

For  10^  the  relations  Vg.  and  8.  with  E!1  «  S*(x)]  =  0 
determine  all  the  jiraps  in  terms  of  any  (non-zero)  one  of 
them. 

The  possibility  of  E]  ^  0  or  of  S' (x)  7^  0  is  covered 
by  the  following: 

The  case  u  +  c  =  0  for  which  10.  reduces  to  dx  =  0, 
corresponds  to  sonic  flov;  which  apparently  is  not  attained 
in  standard  pulse-jets.   The  case  u  +  c  =  F  for  which  10. 
reduces  to  dx  =  Fdt,  corresponds  to  a  sonic  flame  speed  which 
also  apparently  is  not  normally  attained  in  pulse-Jets. 
The  case  u  =  0  7^  F  for  which  9.  becomes  dx  =  0,  introduces 
no  contradictions  even  if  t)'(x)j  7^  0.   The  remaining  special 
case  u  =  F  for  the  characteristic  9.  is  not  exceptional  unless 
E]  7^  0  at  the  flame  boundary  (with  intrinsic  flame  speed 
f  =  0).   In  this  last  case  continuity  of  T  must  be  abandoned 
(cf.  p.  22)  and  the  last  equation  in  7,  must  be  deleted  and 
all  functions  of  t  in  8.  must  be  included  behind  the  jump 
signs  before  using  the  equations  to  determine  the  j\imps  in 
terms  of  a  prescribed  E]  at  the  interface  concerned. 

Having  determined  the  differential  equations  of  the  only 
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curves  (characteristics)  along  which  discontinuities  may 
'propagate*  there  remains  the  problem  of  determining  how 
p,  T  and  u  vary  after  starting  with  prescribed  values  (cf, 
p.  26). 

I.   Some  explicit  analytical  solutions.   For  an  isentropic 
gas  with  E  =  0  _3,  determines  p  as  an  'adiabatic*  f\inction 
of  t  (with  c  =i'/-Pt  ),  and  9^  ceases  to  represent  a  true 
characteristic  since  then  all  the  jumps  across  it  must  be 
zero.   Moreover^   *'    along  the  characteristics 

dx  =  (u  +  c)dt   the  quantities  u  +  o"  respectively  (with 

(T  =  -  /££I)  are  constant  if  the  tube  is  straight  (S(x)const. ). 
If  the  adiabatic  law  in  this  case  takes  the  form  pn  =  const, 
where  Y  =  (2n  +  3)/(2n  +  l)  (with  integer  n  =  o),  then  the 
general  solution  may  be  written  in  the  following  form: 

f  X  -  ut  =  (i  ^)''[G(u  +  a-)  +  H(u  -  0-)] 

11*   K  »  ^ 

1   .  ct  =  <r  (1  ^)''  (-G(u+cr)+H(u-(r)j 

where  G(u+ff)  and  H(u- <r)  are  'arbitrary'  functions  which  are 
determined  by  the  initial  and  boundary  conditions,  and  which 
have  continuous  (n+l)th  order  derivatives  except  for  values 
of  their  arguments  corresponding  to  characteristics  on  which 
actual  discontinuities  arise.  Explicit  determinations  of 
G  and  H  are  possible  if  the  initial  and  boundary  conditions 
are  not  too  involved.  The  method  appears  to  be  most  practical 
if  piecewise  linear  functions  are  used  in  the  conditions  and 
it  may  be  useful  for  example  in  the  discussion  of  the  non- 
linear equations  in  the  region  above  E^  in  Pig,  8,  page  25, 
Details  will  not  be  considered  here. 

II*  General  numerical  solutions.  By  replacing  differentials 
by  differences  (of  sufficiently  small  size),  rather  accurate 
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numerical   solutions  of  the  thermo-hydrodynamical  equations 
1.-5.    subject  to  prescribed  initial  and  boundary  conditions 
may  be  found.     The  type  of  step-by-step  calculation  required 
is  Indicated  by  the  following: 

cjuppose  that  p,  -r  and  u  are  prescribed  at   suitably 
closely  spaced  points  along  a  burve  in  an  (x,t)  diagram, 
with  boundary  conditions  given  at  the  ends.      (Normally  the 
basic  curve  is  the    'initial  curve*  t  =  0  from  which  subsequent 
curve  points  are  found  as  indicated  below.)     To  determine 
points   (and  the  associated  values  of  p,  r  and  u)   on  the 
next  curve  in  the   step-by-step  process,  the  relations  2,-4jj^ 
and  the  differentials 

12.     dp  =  p^dx  +  p^dt   ,  dr    =    T^dx  +      '^^.dt,  du  =  u^dx  +  u^dt, 

may  be  used  thus: 

Let  (x,t)  and  (x  +  6x,  t  +  5t)  be  two  adjacent  points 
on  the  former  curve  and  let  (x  +  a  x  ,t  +  a  t)  be  a  nearby 
point.  In  order  to  keep  track  of  discontinuities  and 
boundaries  it  is  apparently  best  to  choose  the  third  point 
to  lie  at  the  intersection  of  a  discontinuity  line  {a 
characteristic,  a  boundary  or  a  discontinuity  line  for  s'Cx)} 
through  the  first  point  and  a  discontinuity  line  (of  a 
different  family)  through  the  second  point.  It  is  then 
possible  to  calculate  the  values  of  p,  r  and  u  at  the  third 
point  in  terms  of  the  values  at  the  first  two,  using  differences 
in  place  of  differentials  in  12,  (dx  —ax  or  6x,and  dt  ~At  or  6t)# 
The  changes  Ap,  at  and  au  corresponding  to  ax.  At  and 
5x,6t  ,  satisfy  linear  equations  which  can  be  solved  explicitly. 
Details  of  the  method  are  somewhat  lengthy  but  are  not 
difficult  and  so  will  not  be  presented  here# 
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KEFERSWES 

(0)  bee  R.A.E.  (Farnborough,  England)  reports,  especially 
No,  2-1317  s:   "pressure  measurements  on  the  P.Z.G.  76 
flying  "bomb  motor"  by  Wolfe  and  Luck  (1944). 

(1)  See  Wright  Field  Power  plant  memorandum  report 

T  SEPL-5-673-55:   "Preliminary  Testing  of  German  Robot 
Bomb  Engine"  by  Bogert  and  Worth  (18  Nov.,  1944). 

(2)  Aircraft  Engine  Research  Lab.  N.A.C.A.  (Cleveland) 
Mar,  16,  1945  Memorandum  "Performance  Tests  of  Ford 
MX-544  Intermittent  Jet". 

(3)  bee  in  particular  BuAer  (U.S.N.  Engineering  Experiment 
Station,  Annapolis,  Md, )  Jet  Propulsion  Note  No.  20: 
"J.B.-2  flame  motion  data  taken  from  Wright  Field 
high  speed  motion  picture."  By  Edelman  (1  June,  1944). 

(4)  see  Aerojet  progress  reports. 

(5)  General  Electric  Data  Folder  No.  4T  278:  "The 
intermittent  jet  cycle"  by  N.  P.  Bailey  (Sept.  1944). 
Also:  Jet  Propelled  Missiles  Panel  (OSRD  for  Coordinator 
of  Research  and  Development,  U.  S.  Navy  Depstptment) 
Report  No.  92  by  N.  P.  Bailey  and  H.  A.  Wilson. 

(6)  see  eg.  P.  M.  Morse's  "Vibration  and  Sound"  p.  201 
(McGraw  Hill  1936),  and  Rayleigh*s  "Theory  of  Sound" 

§311  and  Appendix  A  (Dover  1945). 

(7)  Applied  Mathematics  Group,  NY  University,  Memorandum 
No.  73:   "Motion  of  gas  in  half -open  pipes"  by 

H.   Lewy  (Oct.   17,   1944). 

(8)  See  eg.  W.   Jost*s  Exploslons-und  Verbrennungs-vorgange 
in  Gasen"  Chapter  III  §  6  (Edwards  Bros.   1943).     Recent 
work  on  flames  and  turbulence  by  Chambre  and  Lin  is 
expected  to  appear  in  the  Joiirnal  of  Aeronautical 
science. 

(9)  see  eg.   N.    P.   Bailey's  "Heat  Engineering"  p.   13 
Wiley  (1942). 

(10)  see  eg.   Lamb's   "Hydrodynamics"  Chap,   I  (Dover  1945). 

(11)  see  NDRC  Applied  Math.   Panel  Report  38. 2R:      "Supersonic 
Plow  and  Shock  Waves"  by  Courant  and  Friedrichs 

(Aug.    1944). 
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(12)  See  reference  (6)  for  "open  end  correction". 

(13)  See  eg.  Dodge  and  Thompson's  "Fluid  Mechanics"  Ch,'  IX 
(especially  p.  213)  (McGraw  Hill  1937). 

(14)  See  eg.  The  linearized  forms  in  Love  and  Pidduck'a 
paper  "On  Lagrange *s  Ballistic  Problem"  T.R, S.  London 
A222,  p.  167-226  (1922),  and  Riemann's  fundamental 
paper  in  his  *Geo.  math,  Werke* ,  p.  145  (Leipzig  1876). 

(15)  See  eg.  P.  de  Haller,  "The  application  of  a  graphic 
method  to  some  dynamic  problems  in  gases",  pp.  6-24 
Sulzei-*  Technical  Review  No,  1  (1945)  (Swiss).  Three 
German  reports,  received  just  as  the  present  AMG-NYU 
article  was  being  printed,  use  analogous  numerical 
and  graphical  methods  (for  isentropic  flow  without 
true  combustion)  .  They  are: 

Schultz-Grunow  ZV7B/PB/RE-2015  -  1  and  2  "Gasdynamische 
Untersuchungen  am  Verpuffungsstrahlrohr".  (1943) 
Sauer  ZWB/PB/RE/1864  -  1  and  2  "Zur  Theorle  des 
Strahlrohrs".  (1943) 

(16)  For  a  general  discussion  of  the  perturbation  series 
method,  see  Frank-Mises,  "Die  Differential  und 
Integralgleichungen  der  Mechanlk  und  Physik,  II", 
p.  127. 


Analogous  numerical  methods  for  cases  involving  flames 
can  be  developed  and  would  be  of  value  for  explicit 
numerical  examples.  Cf.  Appendix  IV. 
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